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We characterize the effect of quantum interference on the line shapes and measured line positions 
in atomic spectra. These effects, which occur when the excited state splittings are of order the excited 
state line widths, represent an overlooked but significant systematic effect. We show that excited 
state interference gives rise to non-Lorenztian line shapes that depend on excitation polarization, 
and we present expressions for the corrected line shapes. We present spectra of 6,7 Li D lines taken at 
multiple excitation laser polarizations and show that failure to account for interference changes the 
inferred line strengths and shifts the line centers by as much as 1 MHz. Using the correct lineshape, 
we determine absolute optical transition frequencies with an uncertainty of <25 kHz and provide 
an improved determination of the difference in mean square nuclear charge radii between 6 Li and 
7 Li. This analysis should be important for a number of high resolution spectral measurements that 
include partially resolvable atomic lines. 



The measurement of accurate atomic transition fre- 
quencies plays an important role in fundamental physics 
from atomic clocks to the determination of nuclear charge 
radii. Determining accurate frequencies requires a suffi- 
cient understanding of the transition line shape. In par- 
ticular, the Lorentzian line shape is of fundamental im- 
portance in the analysis of resonant phenomena in many 
areas of physics fl[. When two or more resonances are 
unresolvable, there arises the possibility of interference. 
The resulting line shape is, in general, no longer a simple 
sum of Lorentzians, even in the low intensity limit. Al- 
though this effect has been known in different contexts 
for many years 043 1 it has typically been ignored in the 
interpretation of Doppler free spectra. We demonstrate 
that quantum interference has an observable effect on 
atomic spectra, which can limit accuracy if not prop- 
erly accounted for. The unresolvable hyperfine struc- 
ture in the D2 lines of hydrogen [(|, lithium [3 potas- 
sium [f|, francium singly- ionized beryllium [l(| and 
magnesium [If[ are examples where interference modi- 
fied Lorentzian line shapes are expected. In this article 
we derive the relevant line shapes, compare them to ex- 
perimental 6 ' 7 Li data, and estimate the systematic errors 
incurred if strictly Lorentzian line shapes are assumed. 



I. DIPOLE SCATTERING LINESHAPE 

We begin with a derivation of the corrected line- 
shape, including quantum interference terms, using the 



* Current address: Department of Physics, College of Science, 
Swansea University, Swansea, SA2 8PP, United Kingdom 

t Current address: Department of Physics, West Virginia Univer- 
sity, Morgantown, WV 26506 

t Current address: Ion Storage Group, Time and Frequency Divi- 
sion NIST, 325 Broadway Boulder, CO 80305 



Kramers- Heisenberg formula [f2[ which describes the dif- 
ferential scattering rate of light incident on an atom ini- 
tially in the state \i) and ending in the state |/). It can 
be derived from Fermi's golden rule [l3| 
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where h is Plank's constant(Zi) divided by 2tt and p s is the 
density of scattered photon states into a solid angle d£l s 
along the scattering direction k s . The scattering matrix 
element Mji is calculated to second order in the electric 
dipole coupling. The scattering matrix element depends 
on the frequency, wavevector and polarization of the in- 
cident light (cjL,kL,£L) and scattered light (cj s ,k s ,e s ). 
The resulting scattering rate is: 
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where c is the speed of light, €q is the permittivity of 
free space, and is the amplitude of the electric field 
of the incident light. The sum is over excited intermedi- 
ate states \j) with transition frequencies u>ji and atomic 
dipole matrix elements Djj = (j|er|z). Here e is the 
electron charge and r is the position operator of the va- 
lence electron. The finite lifetime of the excited states \j) 
are accounted for [f2[ by including the imaginary part 
iTj/2 in the transition frequency Wu, where Tj is the 
inverse lifetime (or full width half maximum for an iso- 
lated Lorentzian line) of Equation [2l valid in the 
low excitation intensity limit, does not include multiple 
scattering effects like optical pumping. Additionally, we 
make the rotating wave approximation, which is appro- 
priate for near resonant excitation. While Eq. [2] is a 
Lorentzian distribution if only one term of the sum is con- 
sidered, since the sum over intermediate states is inside 
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the square, one can see that interference from different 
excited states \ j) is possible. 

For a concrete experimental comparison, we restrict 
our analysis to the case where states \i) and |/) are hy- 
pcrfine states of a single electronic ground state with 
electronic angular momentum J, and the intermediate 
hyperfine states \j) belong to a single excited electronic 
state with angular momentum J'. The states are labeled 
by their total angular momentum and z-projection of an- 
gular momentum \Fi,mi), \Ff,rrif), and \F',m'). 

One can evaluate the atom field coupling matrix el- 
ement by repeatedly applying the Wigner-Eckhart the- 
orem. The reduced matrix elements can be written in 
terms of the electronic excited state linewidth T and a ref- 
erence intensity Iq (see Appendix |A| . (For a closed tran- 
sition such as the Li 2s-2p transitions considered here, 
T = Y j .) This gives 
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ul-w^j?,, and Ap ™ are the normalized 
dipole matrix elements containing all the angular depen- 
dence of the atomic dipole. The explicit form for Af^.™. 
is given in Appendix |A"1 

Since the denominator in Eq. [3] is independent of m', 
we can sum the numerator over to' . Defining the function 

c£/(e s , e L ) = £(e s • Af# )(A#£ • 6 L ), (4) 



we have 
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where Cf_^t (e s , cl) depends on the initial and final state 
quantum numbers Fi, TOj, Ff and rrif. 



Equation [5] describes the differential scattering rate 
of light into solid angle dtt s (along k s ) with polariza- 
tion e s for atoms starting in state \Fi,rrii) and ending 
in | F/, to/). In a typical spectroscopy experiment, the 
final scattering state is unresolved, so the scattering rate 
Rp imi ^F fmf is summed over final states Ff and to/. To 
further simplify the discussion, we assume the detection 
is polarization insensitive and sum over the two scattered 
polarizations e s _L k s for a given detection direction k s . 
If, in addition, we assume an unpolarized atomic sample, 
we must average over all initial . Summing and evalu- 
ating the square in Eq. [5l gives rise to sums of Lorentzian 
components and cross-terms 
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where the line strengths /(k s , cl, F') and cross-term 
strengths g(k s , cl, Fi, F' , F") for a particular laser polar- 
ization and detected direction are given by 



f(k s ,e L: F t ,F') 



s 7 m,i Ffirif 



g{k s ,e lj ,F i ,F',F") = — E C^ f (e s ,e L ) [cC f (e s ,e h ) 

s,rrii FfTUf 



(7) 



where <?t — Si(2-^i + 1) i s the total number of Zeeman 
states in the ground electronic state, assumed here to be 
uniformly thermally populated. When the excited state 
hyperfine splitting is not well resolved, A p. — A p. = 
Ap, w r, then the cross-terms are not necessarily negli- 
gible, as implicitly assumed in the latter portion of [14j . 

Angular dependence: Dipole scattering of light follows 
a dipole radiation pattern (l5| . which for linearly polar- 
ized light depends only on the angle 7 between excita- 
tion laser polarization £l and the fluorescence collection 
direction k s . The angular dependence of the dipole scat- 
tering is proportional to cos 2 7, and it can always be writ- 



ten as a sum of a spherically symmetric component and 
a dipole component (A to t + BP2(cosj))/4:iT. Here A to t 
is the total line strength integrated over all solid angle, 
Pzfa) = (3a^ 2 — l)/2 is the second Legendre polynomial 
(which has zero integral over solid angle), and B charac- 
terizes the amplitude of the angular dependence. By con- 
struction / contains all the scattering linestrength, the 
integral of the cross-terms g, proportional to P2(cos7), 
over solid angle vanishes. A consequence of this angu- 
lar dependence is that /(k s , cl, Fi, F') does not provide 
the correct ratio of line strengths of the Fi — > F' tran- 
sitions for an arbitrary choice of detection direction, 7, 



3 




'z 




J 








y 



FIG. 1. (color online). Coordinate system: The excitation 
laser propagates along x, so that the linear polarization di- 
rection el lies in the y-z plane, parameterized by (9l- The 
detection direction k a lies in the x-z plane and is parameter- 
ized by 9 a . For Doppler-free spectroscopy, the atomic beam 
is taken to be along y. 



width [HGH. Explicit expressions for f^Fi F') 
and g{9\ J ,F i — > F',Fi — > F") are evaluated for lithium 
with the collection along the z direction in Appendix [SJ 
Lineshape Impact on Extracted Frequencies: We now 
give a qualitative discussion of the effect of the addi- 
tional interference terms on Doppler-free, or nearly free, 
spectra. We choose 6 ' 7 Li as an example because of its 
fundamentally unresolvable structure (Ap, /T rj 1) and 
because it allows for direct comparison to experimen- 
tal data. Fig [5] illustrates two primary effects. First, 
the maxima of the total lineshape are shifted relative to 
what is predicted by a simple sum of Lorentzian distribu- 
tions, which can lead to errors in extracting the weighted 
line center. Second, peaks may vary in intensity and 
prominence depending on the polarization angle of the 
laser. For example in Fig[5J #l = 0, the amplitude of the 
F = 2 — > F' = 3 component is reduced with respect to 
the F = 2 — > F 1 = 2 component. 



since B/A to t is not the same for different F 1 . As we will 
show, however, there exist "magic" orientations where / 
does give line strengths consistent with resolved transi- 
tions. More importantly, at these magic conditions the 
cross-terms g vanish, giving rise to purely Lorentzian line- 
shapes. 

We parameterize 7 in terms of angles relevant to an 
experimental geometry. The wave vectors fci, and k s de- 
fine a plane which we take to be the x-z plane. With- 
out loss of generality we can take along x, so that cl 
lies in the y-z plane, making an angle 6*l with respect 
to z, and k s lies in the x-z plane making an angle 9 S 
with respect to z, see Fig. [TJ The scattering is then pa- 
rameterized by the linearly independent angles 9 S and 
0Li /(k s , e L , F,F') = f(9 s , h , F, F') and similarly for g 
jl6| . The spherical harmonic addition theorem [T 7 ] can 
be used to relate P2(cos7) to 9 S and 9^: 

P 2 (cos 7) = - (3 cos 2 6 S cos 2 6 h - 1) . (8) 

The general form for / and g is then 

f(9 s , 6 U F, F') =A F F +^ (3 cos 2 9 S cos 2 9 L - l) 

g(9 s , L , F F', F") = ^— - (3 cos 2 9 S cos 2 9 h - l) , (9) 

where Ap , Bp and Cp F are constants determined 
by evaluating EqJTl When cos 9 S cos 0^ = 1/ v3, g 
vanishes and / correctly gives the line strength ratios. 
This can occur for a range of geometries. In partic- 
ular, when the detection k s is orthogonal to the exci- 
tation kL (6>l = 7, 9 S — 0) as in our apparatus 0, 
then $l = arccos(^) = 9m ~ 54.73° is the so called 
"magic" angle. Similar magic angle effects occur in 
quantum beat spectroscopy, which could be viewed as 
a time domain analogue of the effect considered here, 
where the excitation pulse width replaces the natural 



0.030 


A 0.0170 












0.025 








0.020 


/ V ao145 






0.015 


-1.0 -0.5 0.0 o.5: 


0.010 








0.005 








0.000 









-20 -10 10 20 

ml 

— (MHz) 
In 



FIG. 2. (color online). The scattering rate (or intensity), ^ 
in arbitrary units, of the F = 2 — ^ F' = 1, 2, 3 Doppler-free 
feature in 7 Li with (9l = 0. Red: sum of Lorentzians with 
polarization independent weights, Blue dashed: sum of cross- 
terms, Black: sum of Lorentzians and cross-terms. The laser 
frequency, wl/27t, is the offset from the F — 2 — > F' — 2 peak 
in units of MHz. The inset is the F=2 F = 2 -> F' = 2 peak 
enlarged to show the shift in line center. 

Line centers are typically determined by fitting a sum 
of Lorentzian functions to the observed spectral profile. 
We characterize the effect of cross-terms on line centers 
(both of individual hyperfine components and of centers 
of gravity of composite features) by taking a Doppler- 
free lineshape given by Eq.[7]with cross-terms and fitting 
to it using only Lorentzian functions (amplitude, center, 
offset, linewidth). We then compare the centers given by 
Eq.[7]to the centers extracted from the fit to estimate the 
effect of the cross-terms on measured quantities. From 
Eq. [TJ] (with #l = 7, #s = 0), one can see that the magni- 
tude of the shifts, proportional to the angular dependent 
terms, has maxima at #l = 0, 7r/2 and the sign of the 
effect changes at #l = 9m- This will be experimentally 
verified in the next section. The size of the shifts in Li 
are on the order of 100 kHz to 1 MHz, large enough to 
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completely overshadow effects associated with Doppler 
shifts and optical pumping. 

To provide an estimate for other transitions not explic- 
itly considered here, we imagine atoms with the electronic 
structure of 6 Li or 7 Li with variable hyperhne coupling. 
We consider shifts of individual hyperfine components 
as the hyperhne splitting is varied. We intuitively expect 
that degenerate resonances would not affect the measured 
line position. In the opposite limit, Ap, /r > 1 we also 
expect the line positions to be unperturbed. These two 
limits imply that there must be an intermediate hyper- 
fine splitting that maximally affects the measured line 
positions. We can see in Fig. [3] that this happens when 
Ap, /r is of order one. 
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FIG. 3. (color online). Error in the measured hyperhne split- 
ting (in units of F) for a F = 1/2 -> F' = 1/2,3/2 transi- 
tion as a function of the assumed separation when the the- 
oretically calculated full line shape is fit as the sum of two 
Lorentzians. The dash dotted (solid) curve is for laser po- 
larization Oh — tt/2 (0). The error in the hyperfine splitting 
is greatest where the assumed separation is about 1.3 times 
the natural line width. The red vertical line indicates the ex- 
position of the actual hyperfine splitting for 6 Li. Its vertical 
extent shows the range of errors that can occur when the laser 
polarization is varied between and n/2. 

To get a feel for the apparent shifts of individual 
components as a function of separation, we consider a 
simple analytically solvable lineshape consisting of two 
Lorentzian profiles with splitting A and equal ampli- 
tude. We take line profiles with and without cross terms 
and determine the component positions for each as the 
zero crossings of their first derivatives. We examine the 
difference of the position of the first component in the 
Lorentzian only profile, Xi,, and the position of the cor- 
responding component in the full line profile including 
cross terms, xp, as a function of the splitting A. In the 
limit of distantly spaced resonances, A/r> 1, the dif- 
ference in line centers is xf — xl — r 2 /4A, in agreement 
with the large splitting limit described in Q ■ These shifts 
at large separation have recently been calculated at the 
1 kHz level in meta-stable He [20| and in principle oc- 
cur in muonic hydrogen, although at «100 MHz they are 
much too small to account for the discrepancy between 



proton charge radius values [2lj, 122J . In alkalis with re- 
solvable hyperfine structure, i.e. 87 Rb and 133 Cs, these 
shifts may also appear at the ~10 kHz level which, while 
much smaller than in unresolvable lines is on the order 
of the reported experimental uncertainties [23|, [24| . This 
zero intensity shift may also arise from fine structure in- 
terference, and for Li is ~860 Hz (below our experimen- 
tal uncertainty). These shifts at large separation may 
be particularly insidious because they would only add a 
weak linear dependance to the background without de- 
forming the lineshape as in the case of unresolvable fea- 
tures. 

We also investigate the dependence of an unresolved 
feature's extracted center of gravity on hyperfine sep- 
aration. Using the same procedure, we generate the 
full lineshape, now with three components (F=2 — > 
F'=l,2,3). We vary the splitting via the magnetic dipole 
constant, A 3 / 2 , while fixing the electric quadrupole con- 
stant at the value appropriate for 7 Li. The same qualita- 
tive behavior occurs, producing extracted center of grav- 
ity shifts which are largest when |A 3 / 2 |(oc Ap, ) is of 
order T. However, as seen in Fig. |4j there is now an 
additional feature. Since there are now two resonances 
that can shift relative to a given resonance, the sign of 
the shift can change for a given laser polarization. 
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FIG. 4. (color online). Error in the measured center of gravity 
(in units of T) for a F = 2 — >■ F' = 1, 2, 3 transition as a 
function of the hyperfine A constant when the theoretically 
calculated full line shape is fit as a sum of three Lorentzians. 
The dash dotted (solid) curve is for laser polarization (9l = 
tt/2 (0). The red vertical line indicates the x-position of the 
actual A constant for 7 Li. Its vertical extent shows the range 
of errors that can occur when the laser polarization is varied 
between and 7r/2. 



II. APPLICATION TO EXTRACTING ISOTOPE 
SHIFTS IN LITHIUM 

Having discussed the nature and theoretical implica- 
tions of quantum interference effects on the observed line- 
shape, we apply our theoretical results to experimentally 
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measured spectra taken at multiple laser polarization an- 
gles. Our additional measurement and analysis provides 
a refined determination of the absolute transition fre- 
quencies of the 6 ' 7 Li D2 lines; see Fig. [5] for the com- 
plete 6,7 Li D line level structure. When combined with 
previously measured Dl values Q these new data pro- 
vide an improved measure of the 6,7 Li excited state fine 
structure, 2s-2p isotope shift, and the isotopic difference 
in the 2 P fine-structure splitting, the splitting isotope 
shift (SIS). The SIS provides the best point of compar- 
ison between theory and experiment. We propose that 
the interference effect we describe here is the root cause 
for some disagreements between previous measurements 
in Li [iB-UHl and for the lack of internal consistency of 
the frequency comb based measurement in K Q. 

Improved spectroscopy of the Li D lines is of broad in- 
terest in physics because the isotope shift of these lines 
may serve as a nuclear-model-independent method to 
measure relative nuclear charge radii, which are espe- 
cially interesting in the neutron rich 8 < 9 ' n Li [2t|. Mea- 
sured isotope shifts for the lithium 2s-2p (D lines) [25l — 
[2H |30T | or 2s-3s p3iT[34l ] transitions can be combined 
with precise theoretical calculations [29|, |H, H(| to de- 
termine relative nuclear charge radii of lithium isotopes. 
Additionally, measured D-line transition frequencies are 
used as input for the calculation of species-specific "tune 
in/out" optical lattices for mixtures of quantum degen- 
erate gases [37l - [39| . 
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FIG. 5. (color online). Relevant 6,7 Li level structure. The hy- 
perfine components for the D2 transition have natural widths 
of order the hyperfine splitting. IS (isotope shift) 

Doppler free spectra of the Li D lines were taken at 
different laser polarization angles #l using the appara- 
tus described in 0, H3| and fit using the lineshapes pre- 
sented here convolved with a Gaussian to account for the 
residual Doppler broadening present in the experiment, 



typically 4 MHz. For resolved resonance features with- 
out a polarization dependence, such as the Dl lines, the 
independent fitting parameters are the line center, the 
overall amplitude, a constant background offset, the nat- 
ural width, and the Doppler width. The polarization 
angle of any given data set was fixed. For the unresolved 
fluorescence features, we limited the number of fitting pa- 
rameters by fixing the excited state hyperfine splittings 
to values calculated in [4l| and in agreement with [42j |. 
In addition we fixed the ratio of the unresolved ampli- 
tudes to values given by Eq. [[JJ with numerical values for 
Ap, Bp , and Cp ' F tabulated in Appendix [Bl A small 
correction was made to account for the effect of the finite 
collection angle of the detector (see Appendix [Cj . 

A single fit to five spectra at different laser polariza- 
tion angles demonstrates good overall agreement, includ- 
ing relative line-strengths. Fig. [5] shows the F — 1/2 — > 
F' = 1/2, 3/2 D2 feature of 6 Li (center, wl/2tt «0 MHz) 
and the F = 2 ->■ F' = 1,2 Dl peaks of 7 Li (left and 
right, w L /27r « ±50 MHz). The 7 Li Dl lines enable the 
single fit to multiple data sets because they have no an- 
gular dependence (in general no Dl lines have angular 
dependence), allowing the effect of background light lev- 
els and laser intensity fluctuations to be compensated for 
in multiple spectra taken at different times. The fit to 
these five data sets used only one natural width and one 
(mass scaled) Doppler width. 




FIG. 6. (color online). Amplitude of scattered light, pro- 
portional to Eq. as a function of laser frequency ojl and 
laser polarization angle The laser frequency is offset from 
the 6 Li F = 1/2 ground state by 446799816 MHz (see table U 
for numerical values). The gray scale surface is the complete 
theoretical lineshape including cross-terms and the red points 
are experimental data taken at #l = 0°, 25°, 51°, 75°and 90°. 
The central feature is the F = 1/2 -> F' = 1/2,3/2 transi- 
tions in 6 Li while the two constant amplitude side peaks are 
the F = 1 — > F' = 1, 2 Dl lines of 7 Li. 

To demonstrate the apparent shifts resulting from 
analysis with an incomplete lineshape in measured 7 Li 
D2data, we fit the same spectra taken at different laser 
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polarizations and extract the line centers, with and with- 
out the cross-terms. In Fig. the red points are line 
centers fit without cross-terms and the black points are 
the same data fit with the full theory. The black points 
are self consistent, independent of laser polarization while 
the red points exhibit a strong polarization dependence. 
The fit to the red data is of the form A + £>P2(cos 
The amplitude of the laser polarization dependent shift 
is of order 1 MHz. Near the magic angle 6*l = 9m ~ 54.7 
the Lorentzian fits give the same linecenter as the full 
lineshape. 
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FIG. 7. (color online). Line center of the 7 Li F — 1 — > 
F' — 0, 1, 2 feature fit from experimentally measured spec- 
tra as function of laser polarization angle with respect to the 
collection direction. Frequencies are offset from the 7 Li F=l 
ground state by 446810695 MHz (see tableU}. The black (red) 
data points were extracted by fitting the data to functions 
with (without) interference cross-terms. Error bars represent 
the uncertainties given in tabled] 



A. Discussion of Systematics 

Angular offset: To accurately extract line positions at 
all polarizations, the angle #l between the laser polar- 
ization and the detection optics must be controlled and 
understood. Using a waveplate, we could precisely define 
#l up to a small unknown offset angle Oq. We improved 
our previous estimate of 9q by geometric measure- 
ments made when disassembling the apparatus, finding 
6» = -0.7(10) degrees. 

As a consistency check, we compared the well known 
ground state hyperfine intervals (GHI) to GHI values 
we measure by subtracting optical frequencies at mul- 
tiple angles We note that for small offsets 9q, the 
line shifts near 9^ = 0, n/2 are insensitive to first or- 
der in #o because the derivative of the angular depen- 
dence (oc sin#L cos#l|0 l =7, e 5 =o) vanishes. This is of 
practical utility since data fit at #l = 0,7r/2 with the 
complete lineshape including cross-terms should be accu- 
rate as well as equal to each other. We found that while 
the GHI's derived from measurements of the resolved DI 
lines 0] were consistent with known values (43|, the GHI's 



derived from the unresolved D2 lines at 6*l = 0,7r/2 dif- 
fered from the known values by as much as 30 kHz. This 
disagreement indicates the importance of intensity de- 
pendent shifts on the fitted line shapes when cross-terms 
are significant. 

Intensity dependent shifts: For isolated lines, the fit- 
ted amplitudes are taken to be free parameters and the 
fitted line centers are independent of fitted amplitude. 
As a consequence, the centers of the resolved lines are 
not sensitive to intensity dependent effects like optical 
pumping that modify the line ratios from their theoreti- 
cal values. For unresolvable lines, however, the fitted line 
positions depend on the fixed relative values of / and g 
used in the fit. The unresolvable lines are therefore sensi- 
tive to intensity dependent effects. To explore the impact 
of excitation laser intensity on extracted line centers, we 
measured a subset of spectra at multiple laser powers and 
performed a full optical Bloch equation (OBE) simulation 
of the scattering, including all the ground and excited 
Zeeman levels Q |4|| . We numerically solve the OBE 
with a time dependent Rabi frequency proportional to 
the intensity seen by the atom as it transverses the Gaus- 
sian excitation laser beam. We then numerically generate 
a Doppler free lineshape by calculating the population af- 
ter a Gaussian pulse as a function of laser frequency. At 
the intensities used here and in Q , we find these intensity 
dependent effects are small but important (ss 20 kHz). 
However, we suggest that larger previously reported un- 
certainties (~ 100 kHz) in 39 > 41 K Q ascribed to optical 
pumping could likely be removed by using a lineshape 
that includes crossterms. 

To quantitatively account for intensity dependent light 
shifts and optical pumping effects on the line positions, 
we generate numerical OBE data at several different in- 
tensities and fit the numerical data using the analyti- 
cally calculated line strengths / and g appropriate for 
low intensity. (We confirm that in the low intensity limit, 
the numerical data matches both the expected line po- 
sitions and line strengths.) We then determine the lin- 
ear intensity-dependent line shifts from this numerical 
data, and apply this shift to the measured line positions 
[46j . The laser intensities were determined experimen- 
tally from the relative line strengths of the resolved fea- 
tures taken at different laser intensities. This estimate of 
the intensity is somewhat lower than estimates based on 
measured beam waists and laser power (typically 3.5 mm 
and 3 /iW respectively) but removes uncertainty associ- 
ated with secondary measurements of beam waist and 
power. For most features, the shift was of order a few 
kHz/VW, but for the 7 Li D2 F = 1 -> F' = 0,1,2 
transitions it was as large as 6.7 kHz//iW (for our beam 
waist). The uncertainty in this correction was set equal 
to the value of the applied shift and represents one of 
the largest sources of uncertainty in the experiment. For 
the unresolvable lines considered here, we find that opti- 
cal pumping can have a larger systematic effect than the 
light shifts alone. Future experiments should be careful 
to work at low intensities to avoid these shifts on unre- 
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solvable lines. 

Doppler correction: The correction of the first or- 
der Doppler effect was determined from simultaneously 
recorded forward and reverse beam signals using a corner 
cube to retro-reflect the excitation laser beam. For the 
polarization independent Df lines [7j the systematic con- 
tribution to the uncertainty of this correction is f .4 kHz 
due to imperfections of the corner cube retroreflector. 
Because the retroreflector does not preserve the laser 
polarization, the Doppler correction for the polarization 
sensitive unresolved D2 lines could not be determined us- 
ing the corner cube, and is taken instead from a linear fit 
of correction versus time for resolved components mea- 
sured on the same day. This is necessary because the 
laser alignment drifts slightly over many hours of data 
taking, and results in a larger Doppler uncertainty of 
about f kHz. 

Results: Including the Doppler corrections and the 
power dependent shifts, the GHI values at #l = 0,7r/2 
are in agreement with each other and the known val- 
ues [HI (see Figs. El and EJ)- The value of 9 that min- 
imizes the sin 9^ cos #l angular dependence is consistent 
with the geometrically determined value. The final re- 
ported line positions, shown in table HI represent an aver- 
age over 6*l ■ A representative uncertainty budget is given 
in table HU 



7 LiF=l -> F=2 
9 Q = -3.1" 




9l (deg) 

FIG. 8. (color online). 7 Li ground state hyperfine inter- 
val (F = 1 — >• F — 2) as function of laser polarization angle. 
The measured GHI was determined by subtracting absolute 
measurements of the excited state F = 1 —¥ F' = 0, 1, 2 and 
F = 2 — > F' = 1, 2, 3 features including Doppler and inten- 
sity dependent corrections (see the text). The red (black, 
blue) data is for an angular offset of -3.7° (-0.7° ,2.3°). The 
red (black, blue) curve is of the form Ag a sin(0) cos(#) + GHIo, 
where GHIo is the value measured in [431 ] and Ag is fit to 
the data. The triangular point is data from [?J re-analyzed 
using the procedure described in the present work. Error bars 
represent the uncertainties given in table [J 

Measurements at multiple laser polarizations analyzed 
with the correct lineshape provide an important tool to 
independently estimate systematic errors associated with 
the offset angle 9q. For example, power-dependent shifts 
such as optical pumping can partially cancel the effect of 
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6 Li F=l/2-» F=3/2 : 
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FIG. 9. (color online). 6 Li ground state hyperfine interval 
(F = 1/2 — > F — 3/2) as function of laser polarization angle. 
The measured GHI was determined by subtracting absolute 
measurements of the excited state F = 1/2 — > F = 1/2, 3/2 
and F = 3/2 —yF' = 1/2, 3/2, 5/2 features. The solid 
line is the value measured in [43]. Error bars represent the 
uncertainties given in table [I] 



9q on the lineshape and GHI. Minimizing residuals and 
comparing the GHI near 9 m can still lead to small sys- 
tematic shifts in the line positions. These effects are more 
prominent in 7 Li than 6 Li, and our new determinations 
of the absolute cog transition frequencies differ from our 
previous results [7j, by 83 kHz and 19 kHz, respectively. 
From the absolute frequencies the excited state fine struc- 
ture splitting (Table IIII[) , as well as the 2s-2p IS and the 
SIS (Table [lV|l are calculated and compared to the exist- 
ing literature. As discussed in (3(|, both quantum elec- 
trodynamic and nuclear size corrections largely cancel 
when calculating the SIS. It is, therefore, the most reli- 
able result of theory and has been suggested as a bench- 
mark for testing the internal consistency of experimental 
data. Previously reported results have disagreed with 
each other and with theory far beyond their reported un- 
certainties (Table IDT)) . Our current result resolves these 
discrepancies and is in full agreement with the most re- 
cent theoretical result [4l| . This supports the theory that 
underlies the use of D-line IS's to determine mean square 
nuclear charge radii for short lived Li isotopes. 

Finally, we calculate the difference in the 6 ' 7 Li nuclear 
charge radii using the measured D2 isotope shifts re- 
ported in Table UVl and the Dl shifts reported in 0. This 
serves as a point of comparison amongst different types 
of measurements including elastic electron scattering [47| , 
optical isotope shift measurements on the 3 Si — > 3 Po,i,2 
transition in Li + (48| . and optical isotope shift measure- 
ments of the 2s-3s, Dl, and D2 transitions in neutral 
LifHrJH |334i3 as shown in Fig. [TUJ We calculate the 
difference in nuclear charge radius using Eq. (40) of [2!| , 



5(r c 2 )( 7 ' 6 Li) = <a( 7 Li) - (rt)CU) 



E ) 



C 



(10) 

where (r^)( I Li) is the mean square nuclear charge radius 
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TABLE I. Measured frequencies of hyperfine components and 
centers of gravity (cog) of the 6,7 Li D2 lines0 



7— 


p 


F' 


1_ i. C\J LICIILj y l 1V.L _L _L Zj 1 


L;l L)Z 


Q /9 


K /9 


44o / yyo ( l.UOiyZl) 




9 /9 


■3 /9 
0/ Z 


440 / yyo / o.yuz^zi ) 




9 /9 


1 /9 


44D / yyo / .D / OyZl ) 




1 /9 


*3 /9 
0/ Z 


4^0 / yyouz .1 1 zyio) 


6 Li D2 cog 


1 /9 


1 /9 


4^0 / yyouo.ooo^iD ) 






446799648.870(15) 


7 Li D2 


2 


3 


446809874.895(20) 




2 


2 


446809884.357(20) 




2 


1 


446809890.170(20) 




1 


2 


446810687.873(25) 




1 


1 


446810693.687(25) 




1 





446810696.445(25) 


7 Li D2 cog 






446810183.163(16) 



Unresolved component splittings for the D2 lines are fixed at 
values calculated from the hyperfine A and B constants of [4l| 



TABLE II. Representative uncertainty budget (kHz) 



Uncertainty 6 Li D2 
Component F = 3/2 — !> F' = 5/2, 3/2, 1/2 


Statistical variation 


4 


First order Doppler effect 


10 


Estimate of 9 m 


3 


Laser power dependent shifts^ 


17 


Laser intensity variation 


3 


Hyperfine constant inaccuracy 


2 


Imaging system imperfections 


2 


Magnetic field shift 


< 1 


Reference frequency 


0.089 


Total 


21 



a Optical pumping, multiple excitation recoil, AC Stark shift 

of the i th isotope in fm 2 , E mcas is the measured isotope 
shift in MHz, E = - 10532. 570(- 10532.025) MHz is the 
theoretically calculated isotope shift exclu ding the finite 
size corrections for the D2(D1) transitions [49( and Co = 
-2.4565 MHz/fm 2 @]. The values of the difference in 
mean square nuclear charge radius are —0.708(4) fm 2 for 
the Dl and -0.702(9) fm 2 for the D2 lines. These values 
are self consistent and have the smallest uncertainties yet 
reported. They bring the D-line measurements into full 
agreement with the best values from electron scattering 
and optical IS measurements on 2s-3s and 3 Si — > 3 Po,i,2 
transitions in Li and Li + respectively. 



III. CONCLUSION 

We have reviewed low intensity scattering theory as it 
applies to the spectroscopy of alkali atoms with unresolv- 
able hyperfine structures. We find that the effects of light 
polarization and quantum interference alter the relative 
line strengths and quantitatively affect the extraction of 



e scatt. □ D2 
m Di 

» 2s-3s 



-fr 



-0.2 
-0.4 
-0.6 
-0.8 
-1.0 



[47] [48] [33] [34] [35] [25] [26] [27] [33] [28] [30] [7] This 

Work 

Reference 

FIG. 10. (color online). Measurements of the difference in 
mean square charge radius between 7 Li and 6 Li. The points 
are grouped by type of measurement and are then ordered 
chronologically within different types of measurement. The 
solid black line is the weighted average of the results of refer- 
ences [35l . I47L UK, along with the Dl value from 0] and the 
D2 value from this work. Error bars for the present work rep- 
resent the uncertainties given in table [TV] all other error bars 
represent the uncertainties given in the original references. 



TABLE III. Excited state fine-structure intervals. 
Interval 51 



D Li 2p *P fs 



7 Li 2p 2 P fs 



Splitting 


S (MHz) 


Reference 




10052 


779(17) 


this work 




10052 


799(22) 


Sansonetti 


[7] 


10052 


76(22) 


Brog [50J 




10052 


044(91) 


Walls [27] 




10052 


964(50) 


Noble [28] 




10052 


862(41) 
932(8)£ 


Das [30] 




10050 


Puchalski(theory) 


10053 


310(17) 


this work 




10053 


393(21) 


Sansonetti 


\7] 


10053 


184(58) 


Orth [51] 




10052 


37(11) 


Walls [27J 




10053 


119(58) 


Noble [28] 




10051 


999(41) 


Das [30] 




10051 


477(80 


Puchalski(theory) 



a All Dl values are taken on the same apparatus and reported 
in 

b The uncertainties reported in [4lll represent only the numerical 
uncertainty and do not include any estimate of the size of 
corrections not included in the calculations. 



transition frequencies from data, even in the low inten- 
sity limit. Optical pumping effects at finite excitation 
power can further complicate the lineshape, which we 
account for numerically. This leads to a revised determi- 
nation of the 6,7 Li D2 line frequencies and splitting iso- 
tope shift. We identify several species: H @, 22,23 Na 

39,40,41 K [g[ ; and 221p r |] ; 7,9,ll BeII ^ and 25 Mg jT 



52], 
111, 



for which these complete lineshapes will enable the next 
generation of measurements. 
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TABLE IV. 7,6 Li isotope shifts. 



Transition 


Shift (MHz) 


Reference 


D2 IS 


10534.293(22) 


this work 




10534.357(29) 


Sansonetti [7] 




10533.59(14) 


Walls [2JJ 




10534.194(104) 


Noble [28] 




10533.352(68) 


Das [30] 


SIS^ 


0.531(24) 


this work 




0.594(30) 


Sansonetti [7] 




-0.67(14) 


Walls [2JJ 




0.155(60) 


Noble [28] 




-0.863(79) 


Das [30J 




0.396(9) 


Yan(theory) [36] 




0.5447(1) 


Puchalski(theory) [41] 


a All Dl values 


are taken on the same 


apparatus and reported 


in a 
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Appendix A 



Defining the components of the matrix elements 
(Ap™ ) q for each J — > J' transition 



I a F'm \ 
X-^-Fm )q 



x " 2J ' + 1 {Fm;lq\F'm') x /ff, 



\J2F' + 1 

the dipole matrix elements can be written as 



Fm )1 



V2TT1 [ Frn )q 



(A5) 



(A6) 



Pulling the reduced matrix element (J'||D|| J) out of the 
sum, Eq. [2] can be written in terms of the inverse scat- 
tering rate V and a saturation intensity Iq, 



1 



T = - = 



t 3TTe Q hc 3 (2 J' + 1) 



and 



r 



irhcT 



(A7) 



(A8) 



giving Eq. [31 where u and A are the frequency and wave- 
length of the transition. 



Expressions for the normalized dipole matrix elements: 
The vector components of A are easiest to describe in the 
spherical vector basis A q appropriate for a + , ir and a~ 
light, where 

A x = -(A x + iAy)/V2 

Aq = A z 

A-x = (A x — iA y )/V2. (Al) 

Using the Wigner-Eckart theorem, the dipole matrix el- 
ements are given in terms of reduced matrix elements 
as 

(D F '™') q = ^=== (Fm; lq\F'm'), (A2) 

where (Fm; lq\F'm') is the Clebsch-Gordan coefficient 
for adding \F,m) to \l,q) to get \F',m'). Under the 
assumption that the hyperfine interaction does not mod- 
ify the electronic structure of the state, the .F-reduced 
matrix elements can be written in terms of J-reduced 
elements 

(^||D||F) = (J'||D||J) v / 7f, (A3) 

where the reduced oscillator strength f F for the F-F' 
transition can be written in terms of Wigner 6-j symbols: 

V^T = (-i) F+I+1+J 'V2F + iV2WTi( f f'\}- 

(A4) 



Appendix B 

Calculation of weights f and g: The dipole radiation 
weights / and g are calculated using the expression for 
A Fm ' ( E( l- to determine C F _^ f (Eq. gj), evaluating 
the sums in Eq. [7] and then comparing to the dipole ra- 
diation pattern Eq. [9] Taking k s along z, (i.e. 9 S = 0), 
with the two scattered polarizations e s i = x and e S 2 = y, 
and cl to lie in the z-y plane as in Fig [TJ the terms in 
the sum are given by 

e sl ■ Ap™ = ({Ap™ )i - {Ap™ (Bl) 

e s2 • A F '™' = ±= ((A F F ';f)x + (A F Fr fu) (B2) 

2 A F'm' _ » sm ^L (, K F'm'\ j_(\F'm'\ \ 
e L ■ A Fm - \[ A Fm )l + ( A Fm )-l J 

+ cos0 L (A£>') o . (B3) 

We report line weights and cross-terms for the D2 tran- 
sitions, 2 Si/ 2 — > 2 P3/2i °f alkali atoms and hydrogen 
with 1 = 1/2, 1, and 3/2 in table JvTVTl and EH] respec- 
tively. 

Note that there is no angular dependence to the 
Dl terms, and therefore no dipole dependence (Bp = 

Cp' F " = for Dl). Also note that C F '' F " = C F "' F ', 
physically this is because scattering through F" is in- 
distinguishable from scattering through F' when the F" 
and F' are overlapped within the natural width. 
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TABLE V. D2 weights and cross-terms for 1 = 1/2 applicable 
to H, "Bell 



F" 



Bt 



C 



F',F" 



1/18 
1/12 
5/12 



-1/36 
1/48 
-7/48 



-1/16 



TABLE VI. D2 weights and cross-terms for 1 = 1 applicable 
to 2 H and 6 Li, 28 Na 



F 


F' 


F" 


a f ; 


< 


n F',F" 


1/2 


1/2 


3/2 


8/18 





-4/81 


1/2 


3/2 




10/81 


-1/81 




3/2 


1/2 


3/2 


1/81 





2/405 


3/2 


3/2 


5/2 


8/81 


16/2025 


-14/225 


3/2 


5/2 


1/2 


1/3 


-7/75 


-1/90 



Appendix C 

Collection optics correction: If fluorescence is collected 
over all solid angle there is no polarization dependent 
modification to the lineshape. The equations given in 
the text are valid for light scattered into an infinitesimal 
solid angle. Here we find the modification to the angular 
dependent part of the line weights and cross-terms due 
to the finite numerical aperture of the fluorescence col- 
lection optics. For a given laser polarization 6l, we may 
integrate over the final scattering directions fc s allowed 
by the collection optics (parameterized by # s , </> s ). 

cl = sin 6>l y + cos 6*l z 

- cos 6 S 



i i ' . :j i - . 

k s = sin 6 S cos 4>s x + sin S sin 4> s y + ( 



e L • fc s = cos 7 



(CI) 



Performing the angular integrations over the isotropic 
part, where dfl s — d<j>d cos(# s ), we find 



2ir,e c 



2tt(1 - cos6» c ) = Sq. 



(C2) 



o.o 



TABLE VII. D2 weights and cross-terms for 1 = 3/2 applica- 



ble to 7 ' s 




Li, 21 < 23 ' 34 Na, 


39 > 41 K and 


87 Rb, 7 ' 9 BeII 




F 


F' 


F" 


A F F ' 


B F F ' 


riF'.F" 

F 


1 





1 


1/24 








1 


1 


2 


5/48 


-1/48 


-1/32 


1 


2 





5/48 





-1/48 


2 


1 


2 


1/48 


1/1200 


1/160 


2 


2 


3 


5/48 





-7/120 


2 


3 


1 


7/24 


-7/100 


-7/400 



The angle dependent dipole part is scaled by 



dfl s Bz(cosj) = 7rcos6*csin #c-P2(cos 6>l) 



0,0 



S 2 P 2 (cos L ) 



(C3) 



Here 8c is the half angle of the fluorescence collection 
cone. For determining experimentally relevant fitting 
functions, the ratio of the constant and dipole part is 
important, and we find that the dipole components are 
reduced relative to the constant components as, 



5*2 /So — cos cos 2 



(C4) 



These scaling factors are included as part of the fitting 
functions to account for the numerical aperture of the 
imaging system. Failure to include these scaling factors 
shifts the extracted line centers by « 6 kHz for 9q — 26.6° 
used in this experiment. 
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